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^ Abstract 

^-^ In this report, we present a dynamical systems' approach to study the exact nonhnear wave- 

T-H particle interaction in relativistic regime. We give a particular attention to the effect of wave 

obliquity on the dynamics of the orbits by studying the specific cases of parallel {9 = 0) and 
perpendicular {6 = —txJT) propagations in comparison to the general case of oblique propagation 



Oh 



^ Q =\ — 7r/2, 0[. We found that the fixed points of the system correspond to Landau resonance, and 

that the dynamics can evolve from trapping to surfatron acceleration for propagation angles obeying 



o 






a Hopf bifurcations condition. Cyclotron-resonant particles are also studied by the construction 



£^ of a pseudo-potential structure in the Lorentz factor 7. We derived a condition for which Arnold 



diffusion results in relativistic stochastic acceleration. Hence, two general conclusions are drawn : 
1) The propagation angle can significantly alter the dynamics of the orbits at both Landau and 



these two mechanisms for Landau and cyclotron resonant orbits could become potential candidates 



^ cyclotron-resonances. 2) Considering the short-time scales upon which the particles are accelerated, 

o 

. for problems of particle energization in weakly collisional space and cosmic plasmas. 
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I. INTRODUCTION 

The wave-particle interaction has long been considered a dominant energy-momentum 
exchange mechanism in space and astrophysical plasmas. Beyond the dense and internal 
boundaries of stars and planetary magnetospheres, space and astrophysical plasmas are 
predominantly collisionless, and populated by distributions functions inconsistent with col- 
lisional equilibrium conditions. These plasmas are also believed to be turbulent systems 
described by a conventional energy cascade from large scales to small scales where dissi- 
pation takes place. While fluid theories provide satisfactory descriptions of macroscopic 
quantities at large scales, they are not equipped to explain the plasma physics at smaller ki- 
netic scales, and one needs to include nonlinear kinetic processes (wave-particle interactions 
and wave-wave interactions) for a correct description of these turbulent and collisionless 
plasmas [1]. 

Studies of wave-particle interactions in space and astrophysical turbulent plasmas have com- 
monly fallen under the scope of quasi-linear theory [21 E] • Quasi-linear theory departs from 
linear theory in conserving energy and momentum through the account of the wave-particle 
interaction, resulting in a diffusion process for an ensemble averaged distribution function so- 
lution to a Fokker-Planck equation. However, quasi-linear theory is constrained by a number 
of severe limitations making it inapplicable for plasmas containing large-amplitude quasi- 
monochromatic electromagnetic and/or electrostatic waves. Indeed, the first assumption for 
quasi-linear theory consists in constraining particle orbits to their unperturbed components. 
A second assumption consists in assuming a wave spectrum sufficiently dense so that in- 
terference between modes is destroyed by phase-mixing. Hence, quasi-linear theory is valid 
only when the bandwidth is broad enough to enable resonances to be maintained as particles 
are scattered, and non-linear trapping effects by individual waves are to weak to be taken 
into consideration. 

Due to the inherent difficulties of strong turbulence theories, test-particle methods have 
become a favored tool for the study of wave-particle interaction beyond the constraints of 
quasi-linear theory. Numerous methods have been developed in the context of the radiation 
belts alone, ranging from guiding-center approximation [1], to resonance-averaged Hamilto- 
nian ^, gyroresonance averaged equations [6J and computer simulations taking into account 
approximate and exact dipolar fields [see review [7j and associated references]. However, the 



general case of oblique has often been avoided in favor of the more tractable case of parallel 
of propagation [8J. A strong case can be made for the neglect of oblique propagation and 
nonlinear effects for small amplitude waves since the appearance of a small parallel electric 
field can not trap orbits [9j. However, if the electric field component of the wave becomes 
sufficiently large, such that nonlinear effects can be triggered, then a rich class of orbits can 
result, and the parallel approximation becomes invalid. 

Recent observations of the weakly collisional plasma of the radiation belts are suggesting 
that obliquely propagating waves with Poynting flux two orders of magnitude larger than 
previously observed whistlers waves are commonly generated in the radiation belts and ap- 
pear correlated with relativistic electron microbursts |JOj . These large amplitude waves can 
propagate at large propagation angles {6 < 70°), and possess amplitudes capable of energiz- 
ing electrons on time scales of the order of the milliseconds [11]. If these large-amplitudes 
waves are shown to be a common observational signature in the radiation belts, the conven- 
tional models used to describe the wave-particle interaction using a quasi-linear formalism 
will have to be revisited. It is not only inaccurate to assume that a particle will execute a 
random walk in pitch angle during the course of one bounce period, but as demonstrated be- 
low, a particle can be irreversibly accelerated to relativistic energies in less than one bounce 
period. 

In this report, we investigate the exact nonlinear wave-particle interaction in the relativistic 
regime. The inclusion of relativistic effects is a sine qua non condition for any attempt at 
solving the outstanding problems which have emerged in radiation belts dynamics as well 
as galactic cosmic ray. Our goal is to provide a general framework for the wave-particle 
interaction by using a dynamical systems' approach. Such approach, although lacking the 
level of self-consistency found in numerical simulations, can facilitate the understanding of 
complex systems such as cosmic and space plasmas and therefore provide for an intuitive as 
well as quantitative leap between theoretical models and simulations. 

The report is written as follows. In section 2 we derive the dynamical system as well as its 
fixed points and invariants. In section 3 we treat the special cases of parallel and purely 
perpendicular propagation. In section 4 we study the general case of oblique propagation for 
the cyclotron-resonance case as well as the Landau resonance. Section 5 contains a discus- 
sion of the general framework for the understanding of the wave-particle interaction and the 
effect of oblique propagation in coUisionless plasmas such as the radiation belts. In section 6 



we conclude and discuss studies currently underway to address limitations of the dynamical 
system approach. 

II. DYNAMICAL SYSTEM. 

A. Equation of motion for the general case. 

Our study begins with the equation of motion of a particle in an electromagnetic field, 
as described by the Lorentz equation. The equations of motion can be written as 



dp 

dt 



E(x,t) + -^ xB(x,t) , (1) 



for a particle of charge e, momentum p = m'-)\ and rest mass m. The Lorentz factor, 7, is 
defined in terms of the relativistic momentum as follows : 



V m'^& 

The electromagnetic field configuration consists of a background magnetic field Bq to which 
is superposed an electromagnetic wave given by (5E, (5B) : 

E(x,t) = 5E(x,t) (3) 

B(x.t) = Bo + (5B(x,t) (4) 

The electromagnetic wave vector k is chosen to point in the z direction, obliquely to the 
background magnetic field lying in the y-z plane : 

k ■ Bo = A;fio cos(^) (5) 

(6) 
5B = 6B^± + 6Byy 

with the wave magnetic field components written as 

5Bx = SB sm( kz — ut) 
6By = 6B cos{kz — ut) 

and Faraday's law, expressed in terms of the Fourier components, providing for the compo- 
nents of the electric field 

ckx6E{k,uj)=uj6B{k,uj) (8) 



We can therefore express the dynamical system in terms of the phase velocity f^ = u/k, 
and the variables p$ = m'-jv<^ ; Vti = eSB/mc-f ; Qq = eBQ/mc-f, resulting in the following 
coupled ordinary differential equations : 

Px = Pyilo cos{9) + (j9$ — Pz)^i cos{kz — tut) + Pz^o sin(^) 

Py = -Px^o cos{6) + {pz - p<i.)^i sm{kz - ut) 

[9) 

Pz = —px^o sin(^) + Px^i cos{kz — ut) — PyVti sin(A;2; — ut) 

z = PzV<^/p<i. 

In the classical case, the dynamical system is composed of four equations, the three com- 
ponents of the velocity plus the position coordinate along k. However, in the relativistic 
case the expression for the Lorentz factor must be obeyed and constitutes a constraint on 
the particle's trajectory. We can keep track of this constraint by adding an equation in the 
expression of a dynamical gyrofrequency : 



fio 



d /eBo\ 
dt V mcj J 



-Qr 



pc 



rn^d^ _l_ p2^2 



p 



Vt^Vtipq, 



Px cos{kz — cut) — Py sin(A;2; — cut) j 



(10) 



If we define the constant 6i = Qi/Qq, it is straightforward to see that Qi = 6iQq, and 
similarly, since p$ = ^$(7), the time evolution of this quantity can be written as : 



p$ = — mf$7 



Vtn 



(11) 



In order to simplify the dynamical system, we can eliminate the explicit time dependence of 
the equations by making the following mathematical transformation : 

P'x = Px, Py = Py, Pz = lw{Pz - P<t>), Z = -fu,{z 

for the Lorentz factor : 



'J<pt) 



(12) 



7w I -■ 

^^ 

We can therefore write the equations of motion as follow : 

p'^ = VtQp'y cos{6) - Vtip'^ cos{kz' /-iu,)hw + VtQ{p'J-iyj + p^) sin(6') 
p'y = -9-op'x cos{9) + Vtip'^ sin(A;z77«;)/7«; 

P'zhw = -^oP'x sin(6') + fiip'^ cos(A;2;77^) - nip'y sm{kz' /'j^) - p^ 
z! = p'zV<i>/p<s, 



(13) 



< 



(14) 



If we absorb the Lorentz factor 7^^, into p'^ and k, that is we write p'^ — )■ p'^/'^w and k - 
and write pq> in terms of {p'^,p'y,p'^,z',QQ), we can express the dynamical system as 

p'^ = Vlop' cos{9) - Vtip'^ cos{kz') + no{p'^ + p^) sin(6') 



khui, 



—^oP'x cos{9) + ilip'^ sm{kz') 

-QoP',M0) + ^ii'^)ip'.^os{kz')-p'sm{kz')) 



(15) 



z' 



p',v^/p^ 

where the refractive index is represented as n"^ = c^/f|. The magnitude of the momentum 
is now written as p' = ^p^^ + p',^ + {pjl^)\ hence the Lorentz contraction factor also 
transforms from 7(p) — )■ 7(p'). The dynamical system for the classical case can be recovered 
by setting 7 = 1 and 1/n^ — )> 0. Indeed, the equations are equivalent to the classical case 
under the following transformations : p ^ u and Q — )■ ^2/7 . Hence the main difference lies 
in the time dependence of the Larmor frequencies and the extra term that goes as 1/n^ in 
the p'^ equation. 



B. Representation in terms of (P, a, <J>, z'). 



It is convenient to express the relativistic momentum in spherical coordinates, that is in 
terms of a magnitude p' and phase angles {a, $). This can be achieved by introducing the 
following scalar and vector variables : 



P\\=P- bo 

P± = Po X (P X Po) = P - P||bo 



(16) 



where bo = Bq/Bq. Using these definitions, we can rewrite the momentum p' = {p'^,p'p'^ 
in terms of the pitch angle a and the dynamical gyrophase $, both defined as 



tan(Q;) 
tan(<I>) 



4 

P'n 
P±l 



P'x 



P'±2 P'y cos(^) + P'z sin(^) 
Hence, all three momentum components in the wave frame are written as : 



< 



Px 
P'v 



Pz 



j9'sin(«) cos(<l>) 

p' sin(Q;) sin(<l>) cos(6') — p' cos(a) sin(6') 

p' sin(a) sin($) sin(^) + p' cos(a;) cos(^). 



(17) 



(18) 
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Using the definitions (17) and the representation of the momentum in (18), we can proceed 
.mic 



to write the dynamical system ( 15 ) in terms of the normahzed variables 

^2 



and the function F{a, $, Z), as follow : 
r dP 



(^27 



n 



V. 



2 ' 



kz' ,T = iOt 



dr 



da 

dr 



sin(a) cos($) sm{6)/62 



■ ^ I { sin(a) sin($) sin(^) + cos(a;) cos(6') ]F{a,^,Z) 



n^ 



-—- cos(q;) cos($) sin(6') 
Pd2 

ASsUos^ d) cos(<l> + Z)- sin^ d) cos($ - Z) 
S1S3 



n^ 



cos{9) sin(a) - cos(a;) sin($) sin(6') F(a, $, Z) 



dr V 02-rsm(a) 

+ 7^4^ ( cos' (?) sin($ + Z)+ sin' (^ sin(Z - $) 
tan(Q;) V 2 2 



^i^s cos($) sin(6') 



n^ 



sm a 



Fia,^,Z) 



dZ 

dr 



'52'53-P( sin(Q;) sin($) sin(6') + cos(q;) cos(^) j 



dSs ^ 61626IP 
dr ri^ 



F(a,$,Z) 



e. 



(19) 



(20) 



F(a, $, Z) = sin(a) sin^ (-) cos(<l> - Z) 

+ sin(Q;) cos^ (-) cos($ + Z) + cos(a) sin(6') sin(Z). 

It is easy to observe that we can recover the classical regime by setting ^3 = 7/7^2 = or 
F{a, $, Z) = 0. We now proceed to study some of the properties of the dynamical system. 

C. Fixed Points. 



A common first step in the study of dynamical systems is to find and investigate the 



properties of fixed (stationary) points. The fixed points of the dynamical system (20) are 



defined as the values in (P, a, $, Z, 7), for wliicli (P = « = $ = Z = 7 = 0). It can be 
demonstrated (see Appendix 1) tliat tlie dynamical system, for — 7r/2 < 9 < 0, possesses tfie 
following values for the fixed points : 

P 

2' 

(21) 

1 

7 



7tan(^); 




a^±e±l, 


*^±^; 




Z = 0,it; 






1 


J,- 


2 
C2 


(l + tan^{9)) 



It is already evident from (21 ) that in the case of parallel propagation (^ = 0), the only fixed 
point is that for the trivial case P = 0. The fixed point for the relativistic regime appears 
therefore similar to the classical one for parallel and oblique propagation [12]. The fixed point 
for the relativistic regime will translate into properties found in the non-relativistic regime, 
but also results in different types of structures in their vicinity. For non-zero propagation 
angles, fixed points identify volumes of phase-space composed of physically trapped orbits. 
The trapped orbits could give rise to kinetic distortions in the distribution functions, such 
as beams and temperature anisotropics, as was revealed in the classical non-relativistic case 
[13] . However, because the relativistic equations possess a constraint in the form of the 



Lorentz factor 7, different effects are shown to arise. 



D. Invariants. 



The dynamical system (20 ) also possesses a number of invariants valid for the general case 
of oblique propagation. Knowledge of these invariants is used to construct pseudo-potential 
structures. In turn, these structures provide information on trapped and quasi-trapped 
orbits. 

1. First invariant : Ii. 



Using equations (19) to normalize equation (10), the equation describing the evolution 



of the gyrofrequency can be written as follow 



4^-*.^^ (22) 



for r = kp/moj. Hence, this equation has an exact solution, providing the following constant 
of the motion : 



h = (^3^12 + n2. (23) 

We can write this invariant in terms of the variables P,a,^,Z, and ^3 as follow: 



h = \j5lP^ + 25^Pj52 + 5ln^ (24) 

The conservation of this quantity will indicate the degree to which the constraint for the 
Lorentz factor (pi) is respected in a numerical scheme. 

2. Second invariant : li- 

A second general invariant can be found and expressed in terms of the normalized vari- 
ables as follow : 

/2 = b2[n^ - 1)7 cos(e) - b2P cos(a) + bx sin(^) cos(Z) (25) 

This invariant underlies a fundamental property of oblique propagation. One can indeed 
rewrite the invariant in the form E = irfjc^ ~ P\\, which means that one needs a change 
in the parallel momentum to change the energy. This is a well-known statement resulting 
from the Maxwell-Lorentz invariant quantity E ■ B = 0, since a parallel component of the 
electric field can not be eliminated by any Lorentz translation, while the physics in a frame 
with E\\ = 0, such as in the case of parallel propagation, is no different, therefore, than the 
physics in a frame where E± = E^\ = for which energy is a constant of the motion. 

III. SPECIAL CASES. 

A. Parallel Propagation : ^ = 0. 

The wave-particle interaction problem has overwhelmingly been treated for the special 
case of parallel propagation. Even though we do not present any new result in this section, 
we find it useful to briefly discuss the parallel case as a means of comparison to the general 



oblique case. Setting ^ = in (20), we recover the following dynamical system : 



( dP 6163P , ,^, ^ ^, 



dr 



n^ 



^ = -<5i(53 cos(<l' + Z) + ^ sin(a)F(a, $, Z) 
dr '"^ 



n^ 



d^ 
dr 

dZ 
dr 



-S3 + -^ sm($ + Z) 



62S3P cos{a) 



(26) 



dSs ^ 61626IP 
dr n^ 



F{a,^,Z) 



F{a, $, Z) = sin(a) cos($ + Z). 



In addition to the two invariants Ji and J2, equations (26) also possesses the following 
constant of motion [27] for n'^ ^ 1 : 



{S2P cos{a) - If = 26162 



n'-l 



n^ 



Psin(a)sin(<l> + Z). 



(27) 



Moreover, the existence of physically trapped orbits for 6 = requires that cos(a) = cos($ + 
Z) = 0, hence, a = ^ + Z = 7r/2. However, this conditions results in $ 7^ 0. Aside from the 
trivial case of P = 0, no fixed point exists and the parallel propagation has the particular 
distinction, with respect to oblique propagation, to not possess solutions for which a particle 
could be trapped in Z. 

The parallel case has been studied in both the classical and relativistic regime. The classical 
treatment covered by Matsumoto p^ and Hamza et a/. [12], have shown that one can find 
exact solutions in terms of elliptical integrals. Lutomirski and S'ndan [15] have studied the 
relativistic case showing that similar solutions were also possible. Roberts and Buchsbaum 
[TU] have also treated the relativistic case with a special focus on the case n^ = 1, for 
which a cyclotron-resonant particle was shown to gain energy indefinitely, while for n"^ 7^ 
1, the particle simply becomes phase trapped at cyclotron- resonance with no net gain in 



energy on average. Using the invariant in equation (27) as well as Ji and I2, those results 



can be expressed in terms of a pseudo-potential equation in the parallel component of the 



10 



moinentuiii that we write as y = S2P cos(a) = S2P\\ '■ 



y_ 

2 



-V{y;6i,62,n'^,l2) 



2 
for the function a^y) 



n^-1 



n^ 



2 r 2 



a{yy a{y) 



y^ + n^52 



Myfiv-^ 



(28) 



[n 



^)/{h+y)- Solutions to equation (28) for V(|/; (5i, (^2,''^^, -^2) <0 



are bound states of the system for as far as parallel momentum is concerned. However, this 
only holds for n^ 7^ 1. In the case of n^ = 1, we can easily recover the unlimited acceleration 
found by Roberts and Buchsbaum^^ from the invariants of the motion. Setting ^ = and 



n 



1 for I2, one finds that P\\ is constant. That is, if a particle is at cyclotron-resonance, 



it will remain so forever (or until the wave damps), and gain energy indefinitely. It is 
demonstrated in the remainder of the report, that unlimited acceleration is also possible 
for oblique propagation and that it underlies a specific property of the fixed points. 



B. Perpendicular Propagation : 9 = —tt/2. 

We now investigate the purely perpendicular case, as its treatment will be useful to 
characterize the dynamics for propagation angles that increase towards |7r/2| . The dynamical 



11 



system is written in the following form : 

( dP sin(a)cos($) 5i5-iP 



dr 



da 



So 



+ 



n^ 



sin(a) sin($)F(a, $, Z) 



1 



dr S2P 

6163 



cos(a) cos($) + 6163 sin($) sin(Z) 



+ 



n^ 



cos{a) sin($)F(a, $, Z) 



^* e e r /^N sin($) 

= -53 - 6163 cos z + : ; 

dr 02-rsm(aj 

^1^'^ /^ ^N 5l53COS($) ,^ „N „/ ^ „X 

+ , /, cos(<l> + Z) + \^ . \' cos($ + Z)F{a, $, Z) 



(29) 



tan(a) 



n^ sinfa) 






-(5253Psin(a)sin($) 



d63 61626IP 



dr 



n^ 



F{a,<^,Z) 



F{a, $, Z) = sin(a;) cos($) cos(Z) — cos(a) sin(Z). 



Similarly to the parallel case, the dynamical system (29 ) possesses its own set of invariants 

written as : 

h = 61 cos(Z) + 62P cos{a) 

Jg = 62Pcos{^) sin(a) + 61 sm{Z) + Z + t. 



(30) 



With the fixed point analysis for this particular case showing that no fixed points exists, i.e., 
a particle can not be physically trapped, we make the assumption that the solution for Z 
takes the form of a linear relationship in time : Z = Zq + /3t, with Zq as the initial condition 
and /3 as a constant. Replacing the solution for Z in the invariant I5 results in the following 
expression : 

h = S2P cos(<l>) sin(a) + 61 sin(Zo + f3r) + Zq + /3t + t (31) 

It is therefore evident that for Is = to be true, the term (/3 + l)r must be either zero, or 
compensated by the momentum in x, P^ = Psin(a) cos($), to grow to minus infinity as r 
goes to infinity. In the absence of accessible Landau and cyclotron resonances, the latter 
solution does not appear acceptable. We can qualitatively demonstrate this assumption by 



noting that for r ^ 61, the following approximation must be respected 



tan($) 



S2I3 



T. 
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Hence, either a) 7 — )■ 00, or b) $ — )■ O.In the first case, if 7 — )■ 00, then P — )> 00 as 
well. Hence for /4 to be constant, stationary solutions giving Py = Pcos(a) ~ constant are 
required. Such solutions would necessitate Py ~ 0. Such constraint means that either P^ = 
or Z = 0. But both solutions are unacceptable since they would imply the existence of a 
fixed point, which has been demonstrated to not exist for the special case of perpendicular 
propagation. In the second case, the requirement that $ — ?■ means that since S^P ~ f < c 
is bounded, Z — )■ 0, which is in contradiction with the evidence that Z must be linear in 
time because of a zero parallel electric field. We are therefore left with the assumption that 
/3 ~ — 1, an assumption that can indeed be verified by numerical integration. 
Without any loss of generality, we set Zq = 0, resulting in the solution Z = —r. Using /4 
we find the following solutions for P\\ : 

62Pii=h-6icos{T). (32) 

Similarly, the solution for P^ can be directly found from Is : 

S2P. = h + Sism{T). (33) 

Using those two solutions we can find the exact differential for ^3 : 



"^^^ -^[hcosir) + hsm{r)]dr. (34) 

■'3 '" 



5\ n^ 



Hence, the following solutions for ^3 : 

1 1 25i 



[J5 sin(r) + /4 — /4 cos(r)]. (35) 



51 5m n^ 

We can finally find the exact solution for the last variable in terms of r from the dynamical 
system equation in Z, that is : 



1 9A 

'^2-Pz = \l~p7W\ "^ ^[-^5sin(r) + /4 - /4Cos(r)]. (36) 

We have therefore derived exact solutions for the perpendicular case, based on the exis- 
tence of two invariants and the nonexistence of fixed points. Two limiting cases can be 
deduced from these solutions. If the wave is sustained for long periods, such that the time 
of interaction with the particles Tint ~ l/ew for e ^ 1, the perpendicular propagation 
results in phase trapped orbits with no net gain of energy on average. In the opposite 

13 



case where the interaction would be short-hved such that Tint ~ e/i^, we can calculate 



the average increment in energy during the time of interaction. If we write (35) in terms 
oi E = ■m'jc^, and assume large amplitude, low-frequency waves such that Si/ 62 ~ 1, 
then E/Eq = a/1 + 25^ /(5|n^7Q ~ 1 + v1/c^ and a particle gains energy of the order of 
AE/E ~ v%/c^ for every interaction. Given a prescription in the probability of interac- 
tion P(t>$, At) with an electromagnetic wave of phase-speed f$, one could build a map to 
describe the nonlinear interaction of a particle in a relativistic turbulent plasma composed 
of highly oblique electromagnetic waves. This qualitative analysis for purely perpendicular 
wave applies for particles that do not belong to Landau or cyclotron resonance. 

IV. CYCLOTRON AND LANDAU RESONANCES. 

A. Stochastic acceleration at Cyclotron resonance : w — kuvn = ±5^^0/7. 

The most commonly studied problems of wave-particle interactions have been addressed 
in the context of cyclotron- resonance. However, we demonstrate below that the case of 
cyclotron-resonance contains further intricacies when the general case of oblique propagation 
and nonlinear interaction is treated in the relativistic limit. In order to do so, we construct 
a pseudo-potential function for a particle crossing resonances. 
The resonance condition is written in terms of the normalized variables as : 

7sin^(^) - Pcos(a) cos(e) = ±5/^2. (37) 

Using the resonance condition to replace the expression of Pcos(q;) in /2, we find the fol- 
lowing expression : 

h cos(^) T s = (527(n^ cos2(^) - 1) + 61 cos(Z) sin(^) cos(e). (38) 

If 7 and Z do not have singularities in their derivatives when the resonance condition is 

respected, the following relationship must be satisfied : 

dj _ (5i sin(g) cos(g) .(ry.dZ 

d^ - S^in^cos^-l) '^^^^^ rf7- ^^^^ 

We can find an expression between Z and 7 from the invariant Ji. In order to do so, we 

write the invariant quantity in the following form : 

dZ n^-1 _ P^ + n^ 
~dT 2~ ~ 272 
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n 



272 

77,2 



P \2 P 

- + 1 -2- 
n / n 



272' 



z/ - « 1. 
n 



(40) 



Hence, using equation (40) in addition to (38) we can replace Z and sin(Z) and find a 
pseudo-potential equation in 7 of the form : 



i' 



+ r(7;<5i,52,^) = 



(41) 



for a pseudo-potential written as : 

1^(7; 5l, (52,^) = - 



^2 — 1 n^ 
2 2^ 



X 



5i sin(0)cos(^) 

82 n^ cos2(0) — 1 

-\ 2 



/3i 



A + 1 

1' 



I2 cos(^) =F -5 
52n'^ cos2(6') — (52 

2n 



-7 



/32- /33-7 



(42) 



for the set of constants /3i, /32, /Ss defined as follows : 



^2 



I3i = n^-1 

61 sin(^)cos(^) 

62 n^ cos^iO) — 1 

I2 cos(6') =F s 



(43) 
(44) 

(45) 



52^2 COs2(6') — 62 

If we set the initial conditions Zq = and 70 = 1, we can write /^s = /32 + 1- Taking the 



second derivative of (41), we find the following expression : 



2 ^..2^1 /33/3i(/3i + 1) , 1 (/3i + 1)^ + (/3| - /3|)(/3i + 1) 

7 = ■jPi'y - ■7P2P1 + 



1 

4' 



+ 



+ 



2 72 '2 

3 /33(/3i + l)2 l (/32_^2)(^^^i)2 

2 



-yo 



(46) 



4 y z 7-^ 

This equation can be used to treat the cyclotron-resonance for different limits. We hereafter 
focus on the relativistic low energy case for which 7 = 70+^7, with ^7 ^ 70. Using Newton's 
approximation to express 7"" ~ 7o'"(l ~ ^^^7/70) and setting 70 = 1, we find the following 
forced oscillator equation : 

57 + e257 = A(/3i,/32), (47) 
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for the frequency squared : 



e' = -\fil - /33/3i(/3i + 1) + Im + 1)' + (/33 - /32)(A + 1)] 

+ 3/33(/3i + If - ^[{fii - Pl){P, + Ifl (48) 

and the constant forcing term 

A = -^/3i(/33 - 1) + \h^ + l)(/3i + 3) 

+ \[{^i + 1') + (/3| - PlWi + 1)] - \{Pl - /3|)(/3i + !)'• (49) 

Figure [T] represents the dependence of G as a function of 9 for fixed values of 5i. It is clear 
that for the range of chosen parameters (f$/c ~ .70, ^2 = 1), the oscillations in ^7 can evolve 
from harmonic solutions to hyperbolic solutions as the amplitude of the wave increases. As a 
result of a large wave-amplitude, that is 5i growing, a wide range of propagation angles will 
result in hyperbolic perturbations for a relativistic particle in cyclotron resonance. Figure 
[2] represents the transition from 0^ > to 0^ < 0. As the wave-amplitude increases, the 
particle transit from trapped-orbits to quasi-trapped orbits in phase-space. If the amplitude 
is further increased, the orbit becomes stochastic. Quasi-trapped and stochastic orbits 
are resulting from the wandering of the particle from one cyclotron harmonic to another. 
Hence, the particle gains energy stochastically. This result is an extension of the overlapping 
resonances studied by Smith and Kaufman p!7] for classical regimes. Using a restrictive 
choice of parameters, they found that wave amplitudes of the order of 5i = 6B/B0 > 15 
were necessary to have overlapping resonances. However, our analysis shows that there 
is a window in parameter space belonging to the relativistic regime that allows for the 
overlapping of resonances for amplitudes two orders of magnitude smaller. Similarly to 
the classical case, large-amplitudes translate into a broadening of the phase-trapping cell. 
Trapping cells are also largest for propagations at 9 = 45°. Plotted in figures (IS]) and Q are 
Arnold tongues, that is regions of parameter space (n^, 9, 61, 62) leading to stochastic orbits. 
It is evident from the Arnold tongues that even though the effect described by our analysis 
is purely relativistic, a wide range of parameters can result in stochastic orbits. 
It should be noted that even though the equations presented in this section also apply to the 
case of Landau resonance, the parameter space, in which unstable orbits and overlapping 
can operate, belongs to velocities that must go beyond the speed of light. Therefore, the 
aforementioned result applies specifically to the case of cyclotron-resonances. 
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FIG. 1: Squared frequency B^ as a function of the propagation angle 9 for a relativistic particle in 
cyclotron resonance. Each curves are for different values in the wave-amplitude parameter spanning 
0.01 < 6i < 1. The four bold red lines correspond, from top to bottom, to 5i = (0.05,0.09,0.1,0.3) 

B. Hopf bifurcation at Landau resonance : uj = kuvii. 

A fundamental property of a given dynamical system can be deduced by investigating 
whether the phase-space density and volume is conserved or not. That is, whether or not 
Liouville's theorem applies [I8]. The validity of Liouville's theorem provides the possibility 
to construct distribution functions and follow their evolution in time. Non-conservation of 
phase-space density, either locally or globally, stems from the existence of either attractors 
or non-bounded orbits. Making use of the invariant I2, we compute the divergence of the 
flow in phase-space as follow : 



1 dV 
V dt 


dt 




dP., dPy dP, 

dP.^ dPy dP, 




_7 




7 



+ 



dZ 



(50) 



for the volume in phase-space V and the phase-space vector coordinate ^ = {P^, Py,Pz, Z). 



Hence, equation (50) hints at the existence of an attractor if the volume in phase-space 



shrinks as 7 — t- 00. In the case where the particle's energy oscillates back and forth such 
as for a volume of physically trapped orbits, we can consider Liouville's theorem to apply. 
But it can be shown that such an attractor does exist. A recently published Letter has 
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shown that the attractor arises from a change in parameters that results in the bifurcation 
of the orbits around the fixed points [19J. Indeed, the stabihty analysis in Appendix B 
demonstrates that to every fixed points, combined values in {6,rP), satisfying the condition 
n^ — 1 = tan^(0), correspond to a bifurcation in stability [28j. That is, an orbit, close to the 
fixed point will experience a transition from a (marginally) stable orbit to an unstable orbit. 
We observe that when the condition in parameter space is respected, and for a large enough 
amplitude of the wave magnetic field, the real part of one of the eigenvalues becomes pos- 
itive. This type of bifurcation for pairs of complex conjugate eigenvalues crossing through 
the imaginary axis, is the well-known Hopf bifurcations |20]. 

Represented in Figures ([s]), ([6]), ([t]), ([s]) and ([o]), are typical families of orbits for pa- 
rameters below, equal to, and above the propagation angle at the Hopf bifurcations 
{9c = arctan ^/nP^^^) for a given refractive index n, respectively. The wave parameters 
are chosen for a large-amplitude {6i = 0.1), low-frequency wave (^2 = 0.0696), but similar 
results also apply to frequencies of the order of the gyrofrequency as long as the wave- 
amplitude is sufficiently large to allow physical trapping. We can observe that when 6 < 6c, 
the particle becomes physically and phase-trapped in the phase-space region centered at the 
fixed point. The particle eventually closes unto itself with no net gain on average in energy. 
For 6 = 6c, the particles belonging to the basin of attraction centered around the fixed point 
becomes locked in pitch-angle a and dynamical gyrophase $, and trapped along Z. This 
locking effect results in the divergence of the momentum to infinity under a uniform acceler- 
ation. This mechanism is similar to the surfatron process commonly studied in the physics 
of lasers and in the problem of wave-particle acceleration in astrophysical shocks [2TI - I23] . 
Such an effect is purely relativistic and requires the presence of the Lorentz-invariant paral- 
lel electric field. The violation of Liouville's theorem belongs to volumes composed of these 
surfing and trapped orbits. However, since the surfing acceleration is so efficient, a wave 
would be expected to damp away before considerations for self-consistency and collisions 
are deemed necessary. The case of ^ > ^c manifests itself through the loss of stability of 
the fixed point and the evolution of the attractor into two-dimensional tori. The particle is 
initially trapped in the a, $ and Z plane, but eventually becomes untrapped in Z while its 
orbit never closes. Such a regime in parameter space can as well result in the acceleration 
of particles. Figure [7] shows that despite the incapacity to trap physically the orbits, the 
particle can be accelerated to relativistic levels. It is therefore clear from the above examples 
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that the fixed points manifest themselves differently as a function of the wave obliquity and 
that the propagation angle is a critical parameter for relativistic orbits in the presence of 
large-amplitude waves. 

V. DISCUSSION. 

A. General framework 

A general framework for understanding the wave-particle interaction for a monochro- 
matic wave can be drawn from the previous theoretical analysis of the dynamical system 
presented in this paper. When the propagation is parallel, that is ^ = 0, the electric field 
can be eliminated by making a transformation to the wave-frame, resulting in the particle's 
dynamics being resolved entirely. The particle can be phase-trapped but never physically 
trapped. When the propagation angle increases, the obliquity becomes manifest through 
the appearance of a Lorentz-invariant parallel electric field. This electric field physically 
traps orbits and can result in the creation of a beam parallel to the background magnetic 
field as well as anisotropies in temperature. Indeed, the oblique propagation can provide 
an explanation for kinetic distortions of distribution functions for relativistic energies, in 
a similar manner that it does for the classical case. As the propagation angle increases, 
the stable fixed point {9 < 6c), responsible for trapped orbits and kinetic distortions of 
distribution functions, goes as P = tan(^) and therefore shifts trapped cells to higher 
parallel velocities. If the stable fixed point is too distant from the tail of the distribution, 
no particles will be trapped. This transition from physically trapped to untrapped orbits 
is singularized by the treatment of the purely perpendicular case. For 6 = — 7r/2, as well 
as for particles that do not belong to the basin of the stable fixed point for 6 ^ — vr/2, the 
dynamics of the orbit can be simplified to a back and forth slushing on the wave with no 
net gain in energy on average. 

On the other hand, if the propagation angle reaches the critical value 6 = 9c, at which the 
stability of the fixed point is destroyed by the Hopf bifurcation, the particle belonging to 
the basin of attraction will be accelerated uniformly to relativistic energies. For 6 > 6c, 
a particle initially belonging to the basin of attraction now becomes chaotic and physical 
trapping is lost. 
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In-between the regions of phase-space composed of physically trapped and surfing orbits, 
resides one further source of particle energization. The acceleration in this case originates 
in the cyclotron-resonance and results in stochastic trajectories. The inclusion of obliquity 
as well as the preservation of nonlinearities and relativistic effects, reveal that for a 
given propagation angle, there is a window in parameter space for which a particle can 
be accelerated in a diffusive manner primarily along the perpendicular direction. This 
stochastic process is similar to that of the overlapping resonances for the classical case for 
obliquely propagating electrostatic waves ||17J. The results described in the section above 
consist indeed of an overlapping of resonances, but does operate for wave-amplitudes about 
two orders of magnitude lower than those previously assumed. The explanation for this 
discrepancy with the classical regime is that as the particle gains energy, the dynamical 
gyrofrequency ^o = eBo/m'-fc decreases sufficiently to allow the particle to wander from 
one resonance to another. 

The acceleration mechanisms described above both have the important and interesting 
particularity to operate on short kinetic time scales. The difference is that one operates 
stochastically and energizes particles primarily along the perpendicular direction, while the 
second results in a locking in pitch angle and gyrophase, and accelerates particles coherently 
and primarily along the parallel direction. 



B. Applications to planetary radiation belts. 

The most recent waveforms measured in the radiation belts have revealed an unexpected 
discovery. Large-amplitudes, monochromatic, obliquely propagating, and bursty waveforms 
were not only repeatedly measured in the radiation belts [TDIEIHSS], but appeared correlated 
with electron energization [10] as well as relativistic microbursts events [26]. The correlation 
between chorus waves and electron energization in the radiation belts is not recent, but 
it is suspected that if such waveforms were more commonly present in the radiation belts 
that they could be the dominant trigger responsible for the energization of electrons on 
short time-scales. A study by Yoon[TT] has shown that if one solves the plasma equations 
self-consistently, that such waveforms were indeed capable to accelerate electrons on kinetic 
time scales consistently with the observations. Even though our study lacks the levels of 
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self-consistency provided by the numerical method developed by YoonfTT], we arrive to 
similar conclusions if we choose parameters consistent with the radiation belts measured 
waveforms. If we integrate the dynamical system for a few wave-periods, and with low- 
frequency 82 = 0.1, large amplitude 61 ~ 0.06 and for propagation angles obeying the Hopf 
bifurcations, we find that keV electrons commonly found in the radiation belts could be 
accelerated on the order of the milliseconds to MeV energies. 

However, despite encouraging results, we would like to leave a few words of caution. We 
can not rule that such a mechanism is at play in the radiation belts and the reasons are 
as follows. 1) There is no clear understanding of the origin of the observed large-amplitude 
oblique waveforms in the radiation belts. Before we can pinpoint their origin, it is impossible 
to attempt any self-consistent approach to the current problem. 2) The observations of these 
waveforms are plagued by uncertainties large enough to seriously undermine any attempt to 
determine precisely one or multiple acceleration processes. In the very case of the surfatron 
at Landau resonance, one would need good resolution for the electric and magnetic field 
components of the waves to obtain propagation angles and wave vectors. 3) Finally, the 
wave forms are observed with an electrostatic component and the analysis above needs to be 
conducted with the addition of this compressive electric component. Even though it can be 
shown that the addition of the electrostatic field with the same phase as the electromagnetic 
components of the fields would result in the same condition for the surfatron process, a 
difference in phase would shift the Hopf bifurcation and have non-trivial effects that needs 
to be scrutinized. 

In such context, we can not claim that such a mechanism is at play in the radiation belts, but 
we do suggest that since electrons with keV energies can be accelerated to MeV energies on 
kinetic timescales, that such mechanism could possibly arise in the radiation belts and other 
space and cosmic plasmas who are suspected to be permeated by equivalent large-amplitude 
waveforms 



VI. CONCLUSION. 

We have developed a dynamical system to model the interaction of an ion with an 
obliquely propagating electromagnetic wave in the relativistic limit. We have given a par- 
ticular focus on the effect of the obliquity on the particle dynamics. It was demonstrated 
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that physical trapping of Landau resonant particles could be identified by the fixed points 
analysis. Perhaps the main conclusion of our study is that the wave-particle interaction of a 
single wave demonstrates a rich diversity of mechanisms (acceleration, surfing, stochasticity, 
trapping) for which the propagation angle is an important and critical parameter. Indeed, 
the most telling observation, is that the physics at one propagation angle 9 can be signifi- 
cantly altered for an angle ^ ± e. 

Even though the prime difference between oblique propagations with parallel and perpen- 
dicular propagations, resides in the inclusion of a region of phase-space for which particles 
are physically trapped, we have shown that the relativistic treatment also translates in cou- 
pled values in (^,f<j.) for which particles are accelerated to relativistic energies on kinetic 
time-scales VLqt < 1. Such a mechanism, even though requiring specific wave-properties, 
can be efficient since it operates on short-time scales, and the volume encompassed by the 
attractor is large enough to affect a non-negligible portion of a distribution function. 
Furthermore, it was shown that relativistic effects enhance the cyclotron-resonant stochas- 
tic acceleration. As a result of the overlapping in resonances, particles can wander through 
multiple resonances resulting in a stochastic increase in energy. This relativistic effect is 
of interest, since it provides acceleration for wave-amplitudes lower than those required for 
classical regimes of overlapping cyclotron-resonance. Such mechanism could pertain and 
be more spread than initially assumed in weakly coUisional plasmas where particles can be 
confined for long-time scales. 

It should finally be pointed out that the model we used is not self-consistent, and will 
therefore require corrections in order to take into account the complexity of space and astro- 
physical plasmas. Among these necessary corrections, the departure from a monochromatic 
spectrum to one composed of a bandwidth, appears today as the most fundamental of them 
all. Even though some of the large-amplitude waves recently measured in the radiation belt 
show a significant degree of monochromaticity, the cosmic and space plasmas are mostly tur- 
bulent, and the inclusion of additional waves to confirm or infirm the nature of the processes 
responsible for the acceleration of particles is an inevitable step. However, the dynamical 
system approach offers numerous advantages and the endeavors for greater self-consistency 
can be achieved accordingly. Indeed, the dynamical system for the general case, once fam- 
ilies of solutions have been found, can be used as a background nonlinear solution to the 
wave-particle interaction, upon which corrections, such as addition of waves, changing polar- 

22 



ization, dispersion effects, inhomogeneous background magnetic field, etc., can all be treated 
as perturbations to the family of solutions of the "nonlinear homogeneous" system. Such 
method could be investigated theoretically and numerically, in the similar methodological 
fashion and with comparable tools that Hamiltonian systems were constructed to investigate 
the impact of nonlinear perturbations. 
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Appendix A: Fixed Points Analysis. 

Fixed points of an n-dimensional dynamical system denote stationary solutions for all n 
variables. Fixed points are defined for values of the variables for which the dynamical system 
equations equal zero. In this particular case, fixed points represent orbits of physically 
trapped particles. In order to find the fixed points we proceed as follows. It is clear at first, 
that in order to have 7 = 0, one needs to have F{a, $, Z) = 0. Keeping this in mind, we 



first transform the dynamical system as represented by equation sets (20) into polynomial 
form by using the following change of variables : 

x = e'"; y = e'*; z = e'^; (Al) 

and writing the different trigonometric functions in terms of (x, y, z). For the time evolution 
of P the equation in terms of the {x,y, z) variables result in 

(e^" _ g-i") (gi* + g-**) sin(^) = (A2) 

(x--)(y+-) = (A3) 

X y 

{x'-l){y' + l) = (A4) 
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since x,y,z ^ 0. We apply the same procedure for the remaining three equations of motion 
and we find that for d = the polynomial equation gives : 



sin(^)(x2 + l)(|/2 + i)/52 

-5i5^Px[cOs{e){z^ + l)(y2 + 1) + (^2 _ 1)(^2 _ 1)] = 



(A5) 



the polynomial equation for $ is : 



-4P(x2 - l)yz 



+ sin(^)[2(5iP(a:2 - l){z'^ + l)y - Ax{y^ - l)z/6263] 
+6iP{x' + l)[(y2 + l)(^2 _ 1) ^ cos(^)(i/2 - l)(^2 + 1)] = 



(A6) 



and the polynomial equation for Z resumes as : 



2 cos{9)y{x^ + 1) - sin(^) (x^ - l){y^ 



(A7) 



The equation (A4) has the following solutions : 



X = ±1; y = ±i 



(A8) 



We look at each solution starting with the case x = ±1. Replacing x in equation (A7) 
cancels the second term and results in the following constraint : 



4cos(^)y = 



(A9) 



Since the solutions of this equation are 



cos(^) =0; y = 0. 



(AlO) 



neither of them is acceptable. We are interested in the case of oblique propagation with 
7^ I and by definition y ^ 0. Hence the case x = ±1 does not result in a fixed point for 
the oblique propagation. 



We now look at the second solution which satisfies equation (A4), y = ±i. We replace y 



in equation (A7) and find 



±icos(^)(x^ + l) + sin(e)(x^ 



x^(±i cos(^) + sm{9)) = -{±i cos{9) - sm{9)) 



je 



ie\ 



J9 



x'{±{e'" + e-'') - {e'' - e-'')) = -(±(e*' + e"'') + {e'" - e"*")) 



(All) 

(A12) 
(A13) 
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which results in 



x^ = -e±2*^ 



X = ±ie 



±i9 



(A14) 



To find the value of z we replace the value of y in the equation (A5) which results in the 
cancelation of the first two terms and gives the following result : 



26iPx{z^ - 1) = 



(A15) 



Since by definition x ^ and we are not interested with the trivial case P = the solution 
for this equation is : 

z = ±l (A16) 

The last step consists in finding the value of P for the fixed point, which can be done by 



replacing y = ±i and 2; = ±1 in equation (A7). 



-4P(±i)(±l)(a;2 - 1) + APsm{e)6i{±i){x^ - 1] 
+8 sin(0) (±l)x/(52(53 - 46iP cos(e) (x^ + 1) = 

which can be written as : 



-P = ±2i 



X sin(6') 
x^-1 626^ 



Replacing x with the use of equation (A14), we find : 



6263P = - tan(^) 



(A17) 



(A18) 



(A19) 



Since P is always positive the solution requires ^ < < 0. Transforming back x, y, z to 
the dynamical system variables a, $, Z we find that the fixed points are given by: 



P 



-7tan(^); 



a = ±9± 



TX 



$ = ±- 



TT 



0,7r. 



(A20) 



Using the values in equations (A20) for F(a, $, Z) sets it equal to zero. Hence, equations 



(A20) are the fixed point equations for the dynamical system (20). We can also express the 



fixed points in terms of {px,py,pz,z','y) as demonstrated in table 1. 

Appendix B: Stability Analysis 

The next fundamental step in dynamical system theory is to investigate the equilibrium 
of the fixed points. In order to do so, we apply a basic Lyapunov linear analysis that can be 
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TABLE I: Fixed points in the {px,Py,Pz,Z,j) representation 



PxO 


PyO 


PzO 


Zo 


70 





— p$ tan(6') 


p$ 





1 


yi-f(l+tan2(^)) 





—p^ tan(^) 


p$ 


vr 


1 


yi-^(l+tan2(^)) 





+p$ tan(^) 


P<s> 





1 


^Jl-fil+t^nHe)) 





+P(j) tan(^) 


P* 


vr 


1 


yi-f(l+tan2(^)) 



found in any textbook on dynamical systems. The method is summarized as follows. For a 
dynamical system x = -F(x) possessing a fixed point xq for which -F(xo) = 0, one can make 
a Taylor expansion around the fixed point and keep the first order terms. That is, writing 



dF 



as J^5x 



J5x. We are 



the dynamical for the perturbation 5x and the Jacobian J — ^^ 
left with the task of solving an eigenvalue problem since we can write the solution to the 
linearized equation as 5x = .^jC^'* for the eigenvalues Aj and eigenvectors ^j, assuming the 
eigenvalues are not degenerate. If one eigenvalue Aj > the system is linearly unstable at 
Xo, if not the system is linearly stable at xq. From the expression for 70 in the previous 
appendix, it is clear that a fixed point does not exist for all parameter values of 6 and 
f$. Since 7 > 1, the argument in the denominator square root must obey the condition 
n"^ > 1 + tan^(^). Hence we write the Jacobian for the dynamical system as follows : 

/ Q cos(e) T^i+sm(6>) Q \ 

^270 5270 

cosM 


-!- 



<527o 
- sin(e) _|_ n^^^\_Si 



Si 



^270 



5itan(e)n2-l 



\ 



70 



/ 



where the ± symbols denote the two values of the fixed points in Z = kz'. Solving the 
eigenvalue problem (J — AI).^ = we find a bi-quadratic polynomial function in A that can 
be written as x(A) = A^ + rjiX^ + ^2 = 0, with the constant coefficients rji and ri2 given by 
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the following expressions : 



5i n^-l ,,, cos2(^) 

r/i = -— tan(g) + 

(5270 ^2 5|7^ 

5i / n2 - 1 . sin2(^) sin(^) 



/ n — i ,„, sm y 
^±2sin^ -^^T 



(5|7q v n^ 5i n^ 

r/2 = ^^^sin(^)cos(^) (Bl) 

A close look at the coefficients of equation set (7) shows that all four eigenvalues will cross 
the zero real axis when the condition 

n^ - 1 = tan2(^) (B2) 

is respected. That is, for parameter values corresponding to 7q"^ =0 and resulting in A^ = 0. 



[1] R. Kulsrud, Plasma physics for astrophysics (Princeton University Press, New Jersey, 2005). 

[2] C. F. Kennel and F. Engelmann, Phys. of Fluids 9, 2377 (1966). 

[3] A. Roux and J. Solomon, Ann. Geophys. 26, 279 (1970). 

[4] Y. Omura and H. Matsumoto, J. Geophys. Res. 87, 4435 (1982). 

[5] J. M. Albert, Geophys. Res. Lett. 29, 1275 (2002). 

[6] J. Bortnik, R. M. Thorne, and U. S. Inan, Geophys. Res. Lett. 35, L21102 (2008). 

[7] D. Shklyar and H. Matsumoto, Surveys in Geophysics 30, 55 (2009). 

[8] X. Tao and J. Bortnik, Nonlin. Proc. Geophys. 17, 599 (2010). 

[9] J. Albert, J. Geophys. Res. 115, A00F05, doi:10.1029/2009JA014732 (2010). 

[10] L. B. Wilson, III, C. A. Cattell, P. J. Kellogg, J. R. Wygant, K. Goetz, A. Breneman, and 

K. Kersten, 38, L17107 (2011). 

[11] P. H. Yoon, Geophys. Res. Lett. 38, L12105 (2011). 

[12] A. M. Hamza, K. Meziane, and C. Mazelle, J. Geophys. Res. Ill, A04104 (2006). 

[13] A. Osmane, A. M. Hamza, and K. Meziane, J. Geophys. Res. 115, A05101 (2010). 

[14] H. Matsumoto, Space Sci. Rev. 42, 429 (1985). 

[15] R. F. Lutomirski and R. N. Sudan, Physical Review 147, 156 (1966). 

[16] C. S. Roberts and S. J. Buchsbaum, Physical Review 135, 381 (1964). 

[17] G. R. Smith and A. N. Kaufman, Physics of Fluids 21, 2230 (1978). 

27 



[18] O. Regev, Chaos and Complexity in Astrophysics (Cambridge University Press, 2006). 

[19] A. Osmane and A. M. Hamza, Phys. Plasmas 19 (2012). 

[20] J. Guckenheimer and P. Holmes, Nonlinear oscillations, dynamical systems, and bifurcations 

of vector fields, vol. 42 (1986). 
[21] T. Katsouleas and J. M. Dawson, Phys. Rev. Lett. 51, 392 (1983). 
[22] H. Karimabadi et al.., Phys. of Fluids B 2, 606 (1990). 

[23] A. A. Chernikov, G. Schmidt, and A. I. Neishtadt, Phys. Rev. Lett. 68, 1507 (1992). 
[24] C. Cattell et al.., Geophys. Res. Lett. 35, L01105 (2008). 
[25] P. J. Kellogg et al.., Geophys. Res. Lett. 37, L20106 (2010). 
[26] K. Kersten et al.., Geophys. Res. Lett. 38, 8107 (2011). 

[27] This invariant is the relativistic equivalent found in previous studies, i.e., [15] and |12j . 
[28] The fixed point at Z = vr is unstable. Hence, no physical trapping is possible for 9 < 9c and 

no uniform acceleration can arise for 6 = 9c. 
[29] The stochastic process for the cyclotron resonance case is unlikely to be relevant to the Earth's 

radiation belts since the amplitudes and phase-speed of the waves are much larger than those 

observed, but could be of interest in the context of cosmic rays where growing evidence of a 

structured spectrum suggests multiple acceleration mechanisms. 



28 



1.3 


^^^^^^^^ 


1.2 


/^/^l(^f^i^^^\^^\ 


1.1 


ifii \(. ^^i-JvJ lin\ 


^ 1 


Iv \Vv:~^^^ JIj ) 


0.9 


w\ \\^^^^,^ -^/Mj 


0.8 


■^^Ov^;;-— ^_---<^y 


0.7 


^^^Sis^^^^i^^ 


0.6 


^-iii^^zi^S- 



3 

2.S 

2 

1.5 



1.5 2 

4> 




1 

1 


r^ 





-100 



-50 

z 




FIG. 2: a)Pitch angle a vs dynamical gyrophase ^ for B^ > 0. The particle is phase-trapped. 
b)Pitch angle a vs dynamical gyrophase <^ for 0^ < 0. The particle is quasi-trapped in phase-space 
c) Lorentz factor 7 vs Z for 8\ = (0.05,0.09,0.1,0.3). When 0^ < 0, the orbit is unstable in 7 and 
depart from the forced harmonic oscillation observed for 0^ > 0. d)Resonance condition quantified 
by s(P, 7,a; 0,(52) for the case of 0^ < 0. The particle travels through multiple resonances as it 
gains energy through repeated kicks. 
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FIG. 3: Arnold tongue in the parameter space 



.n 



: c/f#), for di = 0.3, 82 = 1. 
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FIG. 4: Arnold tongue in the parameter space {n,52), for 5i = 0.5, 6 = 45°. 
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FIG. 5: Case 9 < dc = arctan\/n^ — 1. Particle orbit for parameters 6i = 0.1, S2 = 0.0696, n? = 4, 



6r — 1° and initial conditions v' 



xO 



0, v'yo 



-v^tan{9) — 1.6^$, v'^ 



zO 



-v^, Zq = 0. The 



particle is physically and phase trapped. It bounces back and forth in the potential well with no 
net gain in energy. 



32 




FIG. 6: Case 9 = 9c- Three particle orbits seeded with different initial conditions show that the 
attractor is periodic. Parameters 6i = 0.1, 62 = 0.0696, "n? = 4. The orbit is locked in pitch-angle 
a and dynamical gyrophase $, trapped along Z and accelerated uniformly. 
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FIG. 7: Case = 6c- Parameters ^i = 0.1, 82 = 0.0696, n^ = 4. The orbit is locked in pitch-angle 
a and dynamical gyrophase $, trapped along Z and accelerated uniformly. 
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FIG. 8: Case 9 > ^c. Particle orbits for parameters 6i = 0.1, 62 = 0.0696, n? = 9. The attractor is 
lost and can give rise to quasi-trapped orbits in the dynamical phase angles. 
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FIG. 9: Case 6 > ^c-Particle orbits for parameters 5i = 0.1, 82 = 0.0696, v? = 9. Despite the lose 
of the attractor particles can still be accelerated to relativistic energy levels. 
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